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Multiple studies have hypothesized parity is associated with pancreatic cancer risk but obtained conflicting 
results. We conducted a meta-analysis (including a dose-response approach) of current available 
epidemiologic studies to investigate the association between parity and risk of pancreatic cancer. Ten cohort 
studies and ten case-control studies including 8205 cases were eligible for inclusion. The combined RR 
(relative risk) of pancreatic cancer for the parous vs. nulliparous was 0.91 (95% CI, confidence interval — 
0.85-0.97, F — 39.0%, Pjj — 0.01). We observed an inverse association between giving birth to two children 
pancreatic cancer risk with RR of 0.86 (95% CI - 0.80-0.93, F - 8.7%, Ph - 0.36). And no evidence 
supported there was non-linear (P — 0.33) or linear relationship (^"=0.14) between number of parity and 
risk of pancreatic cancer. Findings from this meta-analysis indicate that giving birth to two children has the 
lowest pancreatic cancer risk, mechanism of this protective effect needs further investigation. 



Pancreatic cancer remains an aggressive malignancy that caused the eighth leading cancer- related mortality 
worldwide'. Currently, the only possibility for cure is surgery and only 20% patients have such chance due to 
late detection and diagnosis^. The extremely low 5-year survival rate (<5%) highlighted the urgency to find 
out the etiology of pancreatic cancer. However, in the past few decades, smoking has been confirmed as the 
greatest risk factor for pancreatic cancer and accounts for approximately 25% of cases'''*. The other suspected 
propensities such as alcohol consuming, BMI, physical activities, chronic pancreatitis and dietary factors yield 
disaccord results and need to be further validated^. Notwithstanding the recent genome-wide association studies 
identified several genetic susceptibility loci of pancreatic cancer""'', only explaining —4% of the cases". 

The differences in sex distribution'^ which could not entirely be explained by cigarette smoking and the fact that 
steroid hormone receptors exist in normal and neoplastic human pancreatic tissue'" '"^ plus the findings that anti- 
estrogenic agents inhibit the growth of pancreatic cancer in both animal and human models''' prompted inves- 
tigators to hypothesize that the reproductive factors may contribute to the etiology of pancreatic cancer. Up to 
now, numerous studies''*" have focused on the roles of reproductive factors such as parity, age at first birth, age of 
menarche, age of menopause played in the development of pancreatic cancer, but obtained inconsistent findings. 
Of these factors the parity was the most extensively investigated and yielded conflicting results. The association 
ranged from protective effects'"" to substantial increased risk with the increasing number of birth™. 

The efforts of making a conclusion across these studies hampered by small sample size rendering limited power 
to detect the association, and inconsistent referent groups making non comparability between studies. Therefore, 
we conducted a meta-analysis (including a dose-response approach) of current available epidemiologic studies to 
precisely evaluate the association between parity and risk of pancreatic cancer. And to extensively control the 
misclassification, we divided the number of parity into three groups and combined the corresponding data in each 
group separately. 

Results 

Study characteristics. Using the search strategy described in method paragraph we identified 107 articles, of 
which 80 articles were excluded after review of the title or abstract. 27 full-texts were reviewed and one article was 
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Potentially relevant articles (n=218) 
Date to November, 20 1 3 
Pubmcd(n=103) 
Embasc (n=l 15) 



Duplicates removed (n-1 1 1) 



107 articles further review 



48 article were excluded after title screening: 

Review and editorials (n=9) 

Non-human studies (n-39) 

32 article were excluded after abstract review: 

Not concern about parity and pancreatic cancer risk 



1 ' 



27 articles ftirthcr reviewed 



3 articles were Excluded due to insufficient data 
2 articles were removed for overlapping population 
2 articles were removed for focusing on the number of 
pregnancy not the number of parity 



20 articles finally included 



Figure 1 | Flowchart of study selection procedure of this meta-analysis. 

excluded due to insufficient data^'. Other two studies^^'^^ were drived 
out due to employ the measure of standard incidence rate ratios 
(SIRs) and standard mortality ratios (SMRs) to estimate the risk. 
In addition, two studies^"*'^^ were removed for overlapping popula- 
tion with another longer duration follow-up report^' and larger 
sample size^^. Two studies^*'^' were removed for focusing on the 
number of pregnancy rather than the number of parity. Finally, 10 
cohort studies"-^"'^'-''''-" and 10 case-control studies"'-'=''^'"'-*° were 



eligible for this meta-analysis yielding a combined number of cases of 
8205 (Figure 1). Characteristics of the 20 included studies are shown 
in Table 1. Among these studies, seven were conducted in 
USA'^'"'-"-'^''^-'^'''", two in Canada^"''', seven in Europe"^''''^'^'"'^"'""^", 
one in Japan^', and two in China''^-*', one in Egypt'''. The average of 
following- up duration of cohort studies was 18.65 years, ranged 
from 7.1 to 38 years. Cohort size ranged from 37459 to 1292462, 
and the number of cases varied from 52 to 1959. Almost all studies 
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Table 2 | Summary risk estimates of the association between parous and pancreatic cancer. 
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RR (95%CI) 
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Cohort 
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0.89 (0.84- 


-0.95) 


<0.01 


28.7 


0.10 


Case-control 
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0.95 (0.76- 


-1.18) 


0.62 
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0.02 


Number of cases 
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27 
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0.05 
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0.10 
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Location 
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-0.89) 


<0.01 


16.3 


0.27 


North America 
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0.37 


20.8 


0.18 


*Pvalue for heterogeneity. 

Abbreviations; OR, odds ratio; CI, confidence interval. 
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Figure 2 | Forest plot of parous and pancreatic cancer risk. 
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provided relative risk estimates adjusted for age (n = 19), most 
studies for smoking (n = 14), BMI (n = 8) and few controlled for 
history of diabetes (n = 5) (Table 1). 

Parous vs. Nulliparous. Eleven studies including 38 reports reported 
the results using nuUiparous (number of parity = 0) as reference. The 
summary multivariable-adjusted RR (95% CI) of pancreatic cancer 
associated with parity for parous versus nulliparous was 0.91 (95%CI 
= 0.85-0.97) (Table 2, Figure 2), indicating inverse association 
between parity and pancreatic cancer risk. Moderate between- 
study heterogeneity was observed among studies (P = 0.01, P = 
39.0%). The 38 report-specific RRs of parous versus nulliparous 
ranged from a low of 0.89 (95% CI = 0.84-0.95) after omission of 
report by Teras et al'" to a high of 0.92 (95% CI = 0.86-0.98) after 
omission of report by Steven et al"". The power for the sample size of 
meta-analysis to detect an RR of 0.91 is more than 95%. 

Different number of parity. The results of effects of different 
number of parity on pancreatic cancer risk were presented in 
Table 3. To explore the effects of different number of parity, we 
divided the number of parity into three groups and included the 
studies using number of parity of 0 or 1 as reference. 11 studies 
were assigned to the first group (number of parity = 2), and the 
summary multivariable-adjusted RR (95%CI) of pancreatic cancer 
associated with give birth to two children versus number of parity of 
0 or 1 was 0.86 (95% CI = 0.80-0.93) with P = 8.7% (Ph = 0.36) 
(Figure 3). The second group (>2 and <5) contained 16 studies 
including 28 reports, and our analysis of these reports yielded the 
combined risk estimate of 0.98 (95%CI = 0.89-1.09) with P = 53.7% 
(Ph < 0.01). 9 studies including 11 reports in the third group (S5) 
yielded the combined risk estimate of0.95(95%CI = 0.75-1.20) with 
P = 66.3% (Ph < 0.01). The sensitivity analysis revealed that the 
pooled RRs for the second and third groups were similar before and 
after elimination of individual report. Considering the first group, 
the 1 1 report-specific RRs ranged from a low of 0.82 (95% CI = 0.76- 
0.89) after omission of report by Teras et al"* to a high of 0.92 (95% CI 
= 0.84-1.01) after omission of report by Steven et al"". 

Dose response meta-analysis. The dose-response analysis on the 
number of parity and risk of pancreatic cancer involved 19 studies, 
we did not find a curve linear association between number of parity 
and risk of pancreatic cancer (P = 0.33 for non-linearity, P = 25%, 
Ph = 0.05). 

The combined relative risk of pancreatic cancer for per live birth 
was 0.98 (95% CI= 0.95 to 1.01; P = 0.14 for linear trend) with P = 
66.2% (Ph < 0.01), indicating no linear relationship between number 
of parity and risk of pancreatic cancer (Figure 4). 

Meta-regression. We conducted a meta-regression to 
comprehensively explore the source of heterogeneity of parous 
versus nulliparous. Six factors such as study design, location, 
whether adjusted for smoking, BMI and diabetes, number of cases 
were included in the meta-regression model. In this model, the Adj 
R-squared was 100.00%, and Prob > F was 0.0024, which indicated 
that the model was significant. After 100 times permutation, only 
study location appeared to be significant (adjusted P < 0.01) to 
explain the between-study heterogeneity. 

Subgroup analyses. In the subgroup analysis of parous versus 
nulliparous, significant protective effects of parous on pancreatic 
cancer were observed in cohort studies (RR = 0.89, 95% CI = 
0.84-0.95, P = 28.7%, Ph = 0.10), number of cases of more than 
100 (RR = 0.89, 95% CI = 0.82-0.99, P = 60.2%, Ph = 0.01) and 
European population (RR= 0.82, 95% CI = 0.76-0.89, P = 16.3%, Ph 
= 0.27). In the subgroup analysis of effects of different number of 
parity on pancreatic cancer, the effects still turned out to be null in 
each stratum in both second and third groups. However, significant 
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Figure 3 | Forest plot of giving birth to two children and pancreatic cancer risk. 



4.37 



Study 



RR (95% CI) 



Weight (%) 



MESQUIT 1992 

Cantor 1993 
Ji 1996 
Karlson 1998 

KREIGER 2001 < • 

Skinner 2003 
Silvera 2005 
Teras 2005 
Lin 2006 

LO 2007 — 
Prizment 2007 
Heuch 2008 

Steven 2009 | 
Chang 2010 — ^ 
Zhang 2010 

Lucenteforte 20 1 1 — :*■ 
Duell 2012 

Lee 2013 - 
Overall (l-squared = 66.2%. p = 0.000) 

NOTE: Weights are from random effects analysis 



0.87 (0.72, 
0.89 (0.79, 
1.12 (1 .01, 
1.03 (0.98, 
0.63 (0.46, 
0.91 (0.84, 
0.98 (0.91, 
0.98 (0.96, 
0.98 (0.87, 
0.95 (0.85, 

1.00 (0.91, 
1.12(1.03, 
0.95 (0.93, 
0.86 (0.76, 
1.05 (0.96, 
0.90 (0.80, 

1.01 (0.92, 
0.96 (0.88, 
0.98 (0.95, 



1.05) 
1.01) 
1.24) 
1.09) 
0.86) 
0.99) 
1.06) 
1.01) 
1.11) 
1.06) 
1.10) 
1.21) 
0.98) 
0.98) 
1.14) 
1.01) 
1.11) 
1.05) 
1.01) 



2.19 

4.06 

5.25 

8.35 

0.91 

6.25 

6.64 

10.27 

4.05 

4.63 

5.46 

6.35 

10.24 

3.89 

6.00 

4.34 

5.26 

5.85 

100.00 



1 

.461 1 

Figure 4 | Forest plot of parity number (per 1 live birth) and pancreatic cancer risk. 



2.17 



SCIENTIFIC REPORTS | 4:5313 | DOI: 10.1038/srep05313 



7 



Begg's funnel plot with pseudo 95% confidence limits 




-1 - 



s.e. ot logor 



Figure 5 | Beeg's funnel plot of parous versus nuUiparous. 

protective effects were observed in the subgroups of cohort studies 
(RR = 0.87, 95% CI = 0.77-0.99, P = 45.3%, Ph = 0.10), number of 
cases of more than 100 (RR = 0.87, 95% CI = 0.78-0.97, P = 44.9%, 
Ph = 0.11) and European population (RR = 0.84, 95% CI = 0.76- 
0.93, P = 14.5%, Ph = 0.33) of giving birth to two children (Table 3). 

Publication bias. The result of Egger's test indicated no evidence of 
substantial publication bias for parous versus nuUiparous (P = 0.52), 
and in the first group (P = 0.54), second group (P = 0.68) and third 
group (P = 0.70). And the of Begg's funnel plot of parous versus 
nuUiparous was presented in Figure 5. 

Discussion 

We systematically reviewed ten cohorts and ten case-control studies 
including 8205 cases on the association between number of parity 
and pancreatic cancer risk. The factor became interesting candidate 
exposure related to pancreatic cancer because the difference of gen- 
der distribution and several biologic researches suggested that estro- 
gen had protective effect against pancreatic cancer*' and 
pancreatic cancer might be in case another hormone-dependent 
cancer'""'^. The current meta-analysis indicated that parous was 
inversely associated with pancreatic cancer risk compared with nul- 
liparous. Giving birth to two children has the lowest risk of pancre- 
atic cancer. We did not observe a non-linear or linear relationship 
between number of parity and pancreatic cancer. To the best of our 
knowledge, it is the first meta-analysis deeply mined the effects of 
different number of parity on pancreatic cancer risk. 

The meta-regression analysis revealed that the study location 
might be the major source of between-study heterogeneity. We fur- 
ther carried out subgroup analysis by study design, study location 
and number of cases. And the between-study heterogeneity was lar- 
gely removed when stratified according to study location, which 
further indicated that it mainly contributed to the heterogeneity. 
In the subgroup analysis of parous versus nuUiparous, significant 
protective effects of parous on pancreatic cancer were identified in 
cohort studies, number of cases of more than 100 and European 
population. Similar with the subgroup analysis of parous, the sub- 
group analysis of different number of parity revealed that giving birth 
to two children was inversely associated with pancreatic cancer risk 
in those previously mentioned subgroups. The sample size of cohort 
studies was larger than the case-control studies and most studies 
conducted in Europe also have larger sample size than the studies 
conducted in North- America. So the detection of significant inverse 



associations might be due to dramaticaUy increased statistic power of 
the combined large sample size. And the underlie mechanisms of the 
geographic variation of the effects of parity on pancreatic cancer risk 
is largely unknown and needs further validation and exploration. 

Our meta-analysis indicated parous women had lower risk of 
pancreatic cancer compared with nuUiparous women. Parous 
women are likely to have had longer periods of exposure to high 
levels of circulating estrogens. And animal studies reported that 
the estrogen had the inhibitory effect on the growth of preneoplastic 
pancreatic lesions or transplanted pancreatic carcinoma in rats"'"'', 
indicating that estrogen may have effect against pancreatic cancer. 
However, the discrepancy between giving birth to two chUdren and 
more than two is not fully understood, several potential biologic 
mechanisms might have been proposed. First, there was a progress- 
ive increase in insulin resistance during pregnancy"", and insulin 
resistance was suggested to be independently associated with pan- 
creatic cancer risk*^. Another explanation may be obesity which is 
known as a risk factor for pancreatic cancer*'' and pregnant women 
are more tend to be obese. So the harmful effects of insulin resistance 
and obesity may counteract the protective effect of estrogen against 
pancreatic cancer. 

Our current meta-analysis has several strengths. First, we included 
ten cohorts and ten case-control studies which enabled us significant 
statistical power to detect potential association. Second, to control 
misclassification, we combined different number of parity separately 
rather than used the highest versus the lowest number. Third, we 
applied the model adjusting for most established risk factors there- 
fore largely controlled the confounders. In addition, no publication 
biases were detected, which suggested the entire pooled result maybe 
unbiased. 

Despite these advantages, some limitations of the current study 
should be acknowledged. First we did not have access to the primary 
data from the studies included in this meta-analysis, so could not 
perform additional adjustments for potential important covariates. 
Second, there is a relatively wide range of values for cutoff of highest 
level of number of parity so we could not accurately assign an expo- 
sure value to open-ended category, which might impact our current 
analysis. 

In conclusion, the current meta-analysis indicates that parous is 
inversely associated with pancreatic cancer risk compared with nul- 
liparous, and giving birth to two chUdren has the lowest pancreatic 
cancer risk. The exact mechanism of this protective effect needs 
further investigation. 
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Methods 

Literature search and eligibility criteria. We searched the Pubmed and Embase 
databases from inception to November 2013 targeting the studies investigating the 
relationship between parity and the risk of pancreatic cancer. To undertake a 
comprehensive literature search, we used the combined terms: the MeSH heading and 
keywords relating to pancreatic neoplasms, the MeSH heading "reproductive history" 
and the keywords regarding to the exposure {reproductive factors, parity, birth), 
without language restriction. Reference lists of eligible studies and relevant review 
articles were also scanned to identify further pertinent studies. Studies were 
considered eligible if they met the following criteria: (a) the design of the study is 
prospective cohort study or case-control study; (b) investigate the relationship of 
parity and pancreatic cancer risk; (c) provide or allow the calculation of odds ratio 
(OR), relative risk (RR), or hazard ratio (HR) estimates with 95% CI. We excluded 
non-human, case reports, comparative studies not using an analytical epidemiologic 
design, or studies not reporting analyses of primary data (e.g., letters, editorials, 
narrative reviews) and not providing sufficient data. When multiple studies pertained 
to the same or partially overlapping populations, we used the result with the longest 
follow-up time or the largest sample size. 

Data extraction. For each eligible study, two investigators (BBZ and LZ) 
independently carried out the eligibility evaluation, data abstraction, and quality 
assessment, and disagreements were resolved by consensus. The following data were 
extracted from eligible study: the first author, publication year, geographic region, 
demographics of participants, sample size (number of case and total participants), 
duration of follow-up or study period, covariates adjusted in the multivariable 
analysis, parity number categories, the corresponding OR or RR (with their 95% CI) 
and person-years for each category. For the RR, we chose the one represented the 
greatest degree control of potential confounders. The individual authors were 
contacted with E-mail if the data of interest were not shown in the publications. 

Statistical analysis. Because the absolute risk of pancreatic cancer is low in human, 
the ORs and HRs were considered equivalent to RRs and we simply used RRs 
representing aU of the three measures. To control confounding factors in the greatest 
extent, we extracted the maximally adjusted RRs (95%CI). 

First of all, we evaluated the overall effect of parous compared with nuUiparous, if 
the study considered nuUiparous as reference, we summed up all the parous categories 
(number of parity >0) as ever parity in each study together and treated these different 
categories as different reports. For those studies considering number of parity of 0 or 1 
as reference, we further assessed the effects of different number of parity. We first 
divided the number of parity into three groups (2, >2 and <5, ^5) based on the pre- 
analysis of the data structure across these eligible articles, and then combined the 
corresponding data in each group separately. We pooled the RRs for the overall effects 
of parous and respective effects of different groups in a random effects model 
described by DerSimonian and Laird*^ which takes into account both within- and 
between-study variability. Then we explored the possible non-linear or linear rela- 
tionship between number of parity and risk of pancreatic cancer using a random- 
effects dose-response meta-analysis according to the method described by Greenland 
and coworkers*^'*^. This method requires that the distribution of cases and person- 
years/number of cases and controls, and the RR with 95%CI for at least three 
quantitative exposure categories. For each study, the median or mean level of expo- 
sure category was assigned to each corresponding RR. When the median or mean 
exposure level per category was not presented in the study, we assigned the midpoint 
of the upper and lower boundaries in each category as value of exposure. If the highest 
category was open-ended, we assumed that it had the same amplitude as the closest 
category. Furthermore, a potential cure linear relationship between parity and pan- 
creatic cancer risk was assessed using restricted cubic splines with 3 knots at fixed 
percentiles (10%, 50% and 90%) of the distribution^". A P value for curve linearity or 
non-linearity was calculated by testing the null hypothesis that the coefficient of the 
second spline is equal to zero. 

We quantified the extent of heterogeneity using Q-test^^ and f score^^. We con- 
ducted a meta-regression analysis and subgroup analyses to explore source of het- 
erogeneity and subgroup analyses were performed, if feasible, according to study 
design, geographic region, and number of cases. Sensitivity analysis was also per- 
formed to assess the influence of each individual report on overall estimates by 
sequential removal of individual studies. Funnel plots and the Egger's test^^ were 
applied to examine the publication bias. AU statistical analyses were carried out by 
Stata (version 10.0). The power of our meta-analysis was calculated using PowerV3.0 
(http://www.mds.qmw.ac.uk/statgen/dcurtis/software.html). 
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